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NATIONAL ADVISORY OOMMI^E FOR AERONAUTICS. 

TECHNICAL MEMORANDUM NO. 455. 

NOTE ON RESEARCH WORK BY HELMHOLTZ AND WIEN RELATING 
TO THE FORM OF WAVES PROPAGATED ALONG THE SURFACE 
OF SEPARATION OF TWO LIQUIDS . * 
By J. M. Burgers. 

In 1889 Helmholtz developed, for the calculation of waves 
along the surface of separation of two different liquids, formulas 
which rendered it possible to obtain a little closer approximation 
than that given by the primitive theory. His researches were sub- 
sequently continued by W. Wien, who has given several different 
approximations.** In the dissertations of Helmholtz and Wien, the 
arguments applied, for the most part, to developments in series. 
It is not always easy to see how they arrived at these develop- 
ments and, on the other hand, we do not obtain any idea of the 
degree of approximation of their results. 

Although the importance of these researches has diminished, 

due to the development of strict analytical methods of Levi 

Civita, Struik, Weinstein and others, I wish to take up again 

*"Sur quelques recherches de Helmholtz et de Wien relatives a la 
forme des ondes se propageant a la surface de separation de deux 
liquides." Reprinted from "Rendiconti della R. Accademia Nazion- 
ale dei Lincei," Vol. V, -6th Series, 1st Semester, No. 5, pp. 333- 
338, Rome, March, 1937. 

** Compare H. Von Helmholtz, "Ueber atmospharische Bewegungen" II 
("Zur Theorie von Wind und Wellen" ) , "Sitz. Ber. Preuss. Akad. , " 
1889, p. 761 ( "Wissenschaftliche Abhandlungen, " Vol. Ill, p. 309); 
W. Wien, "Sitz. Ber. Preuss. Akad., 1 ' 1894, p. 509, and "Lehrbuch der 
Hydrodynamik" (Leipzig, 1900), p. 189. 
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the considerations of Helmholtz, "in order to present them in a 
little clearer form which, while not affecting the order of approx 
imation, renders possible a clearer perception of their principles 
and of what they can give. 

Let us adopt a system of coordinates Oxy constantly con- 
nected with the waves, so that the motion of the liquids is sta- 
tionary with respect to this system (Fig. l). Let ABC represent 
the profile of the wave, which we assume to be symmetrical with 
respect to the crests (and hence also with respect to the troughs) 
Let X represent the length of the wave. For brevity, let 
n = 2 tt/\. Let U 1 represent the velocity of the upper liquid 
in the direction of the positive x's, and U 3 the velocity of 
the lower liquid in the opposite direction. 

Let z = x + iy. Let x x = ^ + i^ x , the complex potential 
of the motion of the upper liquid, and x. s = cp 2 + 1% , the simi- 
lar function for the lower liquid. The first of these two func- 
tions gives us a conformal representation of the part situated 
above the profile of the wave ABO of a band of the width X of 
the plane of the z's on a semi-band of the width T3 1 X of the 
plane of the function x j( The second function gives us the con- 
formal representation of the lower part of the same band of the 
plane of the z's on a semi-band of the width XL, X of the. plane 
x 2 (Compare Figs. 2a and 2b. I n all the figures, the correspond- 
ing domains are designated by similar hatchings). The immediate 
determination of these two correspondences presents difficulties 
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due to the unknown shape of the profile ABC. The methods of 
Helmholts and Vv'ien consist in replacing the true profile by anoth- 
er curve of such shape that it is possible to make it correspond 
to a straight line "by means of a simple auxiliary function. 

The function introduced by the above-mentioned authors through 
the formula 

„ , . -inz ny-inx ,., v 

v=a+iT=e = e J (1) 

gives a transformation of a band of the plane of the z's in the 
whole domain of the plane of the variable v. The profile of the 
wave is transformed into a closed curve (Fig. 3), For y = + 
v becomes equal to °° . For y = - & , v becomes equal to zero 
(the point F in Fig. 3). In the case of waves of infinitely 
small amplitude, the curve ABO in the plane of the v' s becomes 
a circle. In order to arrive at a more general case, Helmholtz 
substituted an ellipse for the true curve. The degree of approxi- 
mation obtainable in this way depends on two parameters: the ec- 
centricity of the ellipse and the location of the center. (The 
dimensions of the ellipse are of no importance, because a homo- 
graphic expansion with reference to the center of the plane of 
the v' s only produces a displacement of the plane of the z's 
parallel to the axis of the y's.) 

In ...agreement, with Helmholtz and 1/Vien, let us put 

v = cos t - E (2) 
in which £ = I + i t] , while E is a real constant, which we 
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can represent by cos c. Hence the equation rj = constant = h 
produces an ellipse in the plane of the v's. The coordinates of 
the points A and B determine the height of the waves in the 
plane of the z's and we readily find . ■ 

v v ~ 2l i o- coshh + cos g 
*B ^A ~ 2 tt ^ coshh - cos e * 

Formula (2) gives a conformal representation of the plane of 
the' v's and hence of a band of the plane of the z's on a semi- 
band of the plane of the £' s (Fig. 4). On recalling the rela- 
tions existing between z and the functions 'x 1 and x 2 , we find 
that the part of the semi-band in the plane of the £' s situated 
above the straight line ABO corresponds to a semi-band of the 
plane of the x^b, while the rectangular part A'B'C'CBAA' of 
the plane of the £ corresponds to a semi-band of the plane of 
the x 2 ' s. 

It is easily found that the first of these two correspond- 
ences is given by the homographic transformation 

x x = t, + constant. (3) 

The correspondence between the domain A'B'C'CBAA' of the 
plane of the £' s and the semi-band of the plane of the x 2 's 
can be formulated by means of an auxiliary variable t, by em- 
ploying the method of Schwarz-Christof f el. As Levi Oivita has 
kindly shown me, it is advantageous to consider only half of each 
of these two figures, obtained by dividing each one by a vertic- 
al line passing through the points B and B' . Let us take the 



H.A.C.A. Technical Memorandum No. 455 5 

left half A' B ! BAA 1 of the domain in the plane of the £»s, cor- 
responding to the right half of the semi-band in the plane of the 
ota's. Let us make these two parts correspond to a semi-plane of 
the variable t (Fig. 5), so that the point A' corresponds to 
t = 0, B 1 to t = 1, B to t = k~ 8 and A to t = «. On taking 
k as the modulus of an elliptical function, and on designating by 
K and K 1 the complete integrals of the first form, the corre- 
sponding parameter being q = e~ , we have 



t = sn* (4) 



On the other hand, if t = c 2 represents the point which 

must correspond to the point F(ty s = + cu ) of the plane of the 
"x. s 1 s , we find 

d-v const. 

2 _ 



dt (t - c 2 ) Jt - k- B 
whence, on suitably determining the multiplicative constant 

t + c s - 2k~ 2 

C0S uT = 1 - o- (5) 

The relation between c 3 and E = cose is given by the 
equation 

c = sn ( . 
V 71 / 

In the figures, it is assumed that E < 1. Then e is real. 

■ The relations between the variable z and the potentials 
Xj and x 2 , which are given by the developments in series by 
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Helmholtz and Wien, now find themselves expressed in finite form. 
The equation which expresses the equality of the pressures 



(P 2 - P x )gy + £ P 2 



may 


be 


written 




2 


1 


d x 


dz 




a p i 


dz 



= constant (6) 



Let us abbreviate by making 



nx = X, ny = Y, nz = Z; P = 



2(P 2 - P x )g 



Q = 



np 2 U 2 2 



2(P 2 - Pjg 



Then, on introducing the variable £, equation (6) takes the 
form 

The derivatives , are real on the line ABC. We 

dC d£ 

can make all the necessary differentiations and, after eliminat- 
ing t, obtain 



| lg |(cos t coshh - E) 2 + sin 2 £ siiiih| + 0 



(7) 



fl~c 2 k 2 sn* TT 



0 



k ( cost coshh-E) + sin 3 t sinh 2 h 
sin i cosh 3 h+cos s £ sinh 2 h 



Equation (7) contains five parameters: k, c, P, Q, and 0. 
It is evident that we cannot satisfy this equation identically, 
i.e., at all the points of the line ABC. We can only demand 
that it be* satisfied in a limited number of points (five at most). 
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In order to retain the possibility of considering different types 
of waves, this number can be further diminished to four. Then one 
of the parameters remains arbitrary. For these points, I have 
taken those which correspond to the values of i, namely,, o, 
tt/3, 2tt/3 and tt . For a series of values of k, or what 
amounts to the same thing, of the parameter q, I have made 
numerical calculations, whose results are recorded in the follow- 
ing table: 



q 


E - cose 


P 


Q 


Ymax — Yjnin 
2tt 


A 


0 


0.828 


0.0678 


0.4322 


0 




0.207 


0.816 


0.0657 


0.480 


0.107 


0.008 


0.307 


0.807 


0.0313 


0.547 


0.156 


0.021 


0.403 


0.800 


0.0528 


0.657 


0.199 


0.050 


0.490 


0.795 


0.0418 


0.831 


0.235 


0.088 


(0.600) 


(0.794) 






(0.279) 





The A column gives the difference between the values of the two 
members of equation (7) at the point i = tt/2, expressed in frac- 
tions of the value of "? max ~ "^rtiin' Tlie Profiles of the corre- 
sponding waves are represented in Figure 6, in which we have 
marked the points where i has one of the values o, tt/3, tt/2-, 
2tt/3, tt . 

As soon as cos e cos h becomes smaller than unity, which 
occurs for q a little larger than 0.49, the profile of the wave,' 
given by formulas (l) and (2), presents a secondary maximum in the 
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troughs, accompanied "by a minimum .on each side (Compare the bottom 
curve in Figure 6, in which the minimum has also been marked). In 
this case the approximation method seems to be no longer practica- 
ble. • 

For small values of q we can develop in series the magni- 
tudes which figure in equation (7). These series can be arranged 
according to the cosines of the multiples of i . If we retain 
only the terms of the order q 3 , the series terminate at cos 3£ 
and we can demand, as was done by Helmholtz, that the coefficients 
of cos o i, cos I, cos 2 i and cos 3 I be equal in both mem- 
bers of equation (7). In this way we again obtain four equations 
between the five parameters, which renders it possible to express 
them, in terms of one of them alone. It was In this way that the 
values of E, P and Q for q = o were found. 1 found that the 
equations thus obtained were the same as the ones given by Helm- 
holtz in his memoir (excepting for one of the coefficients which 
figure in the development of Y ). 

/ 

Translation by Dwight M. Miner, 
National Advisory Committee 
for Aeronautics. 
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Fig. 2a Fig. 2b 
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Figs .4, 5&6 





q=0.307 

q=0.307 4 
q=0 . 403 i>- 
q=0 . i90 ^ 
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x=0 




Fig. 6 



